A frustrated quantum spin-s model on the Union Jack lattice 
with spins s > \ 



R.F. Bishop and P. H. Y. Li 

. School of Physics and Astronomy, The University of Manchester, Schuster Building, Manchester, M13 9PL, United Kingdom 
Received: date / Revised version: date 

Abstract. The zero-temperature phase diagrams of a two-dimensional frustrated quantum antiferromag- 
netic system, namely the Union Jack model, are studied using the coupled cluster method (CCM) for the 
two cases when the lattice spins have spin quantum number s = 1 and s = | . The system is defined on a 
square lattice and the spins interact via isotropic Heisenberg interactions such that all nearest-neighbour 
(NN) exchange bonds are present with identical strength Ji > 0, and only half of the next-nearest-neighbour 
(NNN) exchange bonds are present with identical strength J2 = kJi > 0. The bonds are arranged such that 
on the 2x2 unit cell they form the pattern of the Union Jack flag. Clearly, the NN bonds by themselves 
(viz., with J2 = 0) produce an antiferromagnetic Neel-ordered phase, but as the relative strength k of the 
frustrating NNN bonds is increased a phase transition occurs in the classical case (s — > 00) at = 0.5 
to a canted ferrimagnetic phase. In the quantum cases considered here we also find strong evidence for a 
corresponding phase transition between a Neel-ordered phase and a quantum canted ferrimagnetic phase 



at a critical coupling k c1 = 0.580 ± 0.015 for s — 1 and « C1 = 0.545 ± 0.015 for s 



In both cases the 



ground-state energy E and its first derivative dE/dn seem continuous, thus providing a typical scenario 
of a second-order phase transition at k = k c1 , although the order parameter for the transition (viz., the 
average ground-state on-site magnetization) does not go to zero there on either side of the transition. 
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75.30.Kz Magnetic phase boundaries - 75.50.Ee Antiferromag- 



> ' 1 Introduction 



In a recent paper [T] we have used the coupled cluster 
method (CCM) [21120] to study the magnetic order on 
a frustrated Heisenberg antiferromagnetic system defined 
on the Union Jack lattice described in Sec. O below. In 
the earlier paper [T] we studied the case of particles with 
spin quantum number s = h. In the present paper we 
further the investigation of this frustrated Union Jack 
model by replacing the spin-^ particles with particles with 
higher values of s. In particular we study each of the cases 
s = 1 and s = | , both of which are computationally more 
challenging than the previous case with s = ^. Just as 
in the earlier paper, however, we again use the much- 
studied CCM. Our main rationale for the present work 
is that one knows in broad terms that the value of the 
spin quantum number s can play an important and often 
highly non-trivial role in the behaviour of strongly cor- 
related magnetic-lattice systems. They often exhibit rich 
and interesting phase structures that are themselves a con- 
sequence of the interplay between the quantum fluctua- 
tions and the frustration due to the competing interac- 
tions present in the system under study. The strength of 
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the quantum fluctuations can itself be tuned, for exam- 
ple, either by introducing anisotropy terms in the Hamil- 
tonian (and see Refs. [SJEMZ] for examples) or by varying 
the value of the spin quantum number s of the particles 
(and sec Ref. [5] for an example). In the present paper we 
study the effect on the Union Jack model of increasing s, 
in order to throw more light on the mechanisms for mag- 
netic ordering inherent in the system. Alternative ways of 
varying the quantum fluctuations, not studied here, would 
be to introduce anisotropy into one or both of the mag- 
netic Heisenberg bonds present, either in real space (and 
see Ref. [7] for an example) or in spin space (and see Refs 
[BUS] for example). 

The general field of quantum magnetism at zero tem- 
perature, for spins on the sites of regular lattices in two 
spatial dimensions [51 ll0l[TTj . has become an important 
and fascinating subject in recent years, and one that is 
at the forefront of modern condensed matter research. 
Much attention has been focussed on frustrated systems 
where different types of bonds are in competition with 
one another, especially when each type acting alone pro- 
duces a different ground-state (gs) ordering and a different 
phase. A set of touchstone problems in this respect is the 
case of the two-dimensional (2D) square-lattice Heisen- 
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berg antifcrromagnet (HAF) with nearest-neighbour (NN) 
exchange bonds with strength J\ > that by themselves 
produce a Neel-ordered phase, but now with the addi- 
tion of frustrating next-nearest-neighbour (NNN) bonds 
of strength J2 > on some or all of the diagonals of some 
or all of the fundamental square-lattice plaquettes. A re- 
view of the "pure" system for J2 = is given in Rcf . [12] . 
Although quantum fluctuations certainly act to destroy 
the perfect Neel antiferromagnetic ordering of the classi- 
cal model (equivalent to the limiting case s 00), it is 
well established that for the s = \ case the order parame- 
ter, namely the sublattice or staggered magnetization, has 
a value equal to about 61% of the classical limiting value. 
For such unfrustratcd 2D models the most accurate results 
are generally provided by quantum Monte Carlo (QMC) 
simulations (and see, e.g., Ref. [T5] for the spin- half HAF 
on the 2D square lattice). 

Once NNN bonds are added to the above pure spin- 
i HAF the situation becomes much more complicated. 
The prototypical model in this respect is the so-called J\- 
J2 model in which all possible NNN bonds are included, 
and which finds good experimental realization in such lay- 
ered materials as Li 2 VOSi04, BaCdVO(P04)2, and oth- 
ers. Various approximate methods have been used to sim- 
ulate the properties of this system, including the cou- 
pled cluster method (CCM) 0[71[T1[IMIS] ' series ex P an " 
sion (SE) techniques [TTllTMTil20"ll2T] . exact diagonaliza- 
tion (ED) methods [221I231I24] . and hierarchical mean-field 
(MF) calculations [23] ■ For frustrated spin-lattice models 
in two dimensions both the QMC and ED techniques face 
formidable difficulties. These arise in the former case due 
to the "minus-sign problem" present for frustrated sys- 
tems when the nodal structure of the gs wave function 
is unknown, and in the latter case due to the practical 
restriction to relatively small lattices imposed by com- 
putational limits. The latter problem is exacerbated for 
incommensurate phases, and is compounded due to the 
large (and essentially uncontrolled) variation of the results 
with respect to the different possible shapes of clusters of 
a given size. 

Several other models of this same general class of spin- 
i models on the 2D square lattice with both NN and NNN 
interactions present, but in which some of the NNN J 2 
bonds are removed, have prompted considerable recent in- 
terest. One such is the Shastry-Sunderland model {261I2T1 
[28] . in which three-quarters of the J2 bonds of the J1-J2 
model are removed in a particular arrangement so that 
each lattice site is connected by four NN J\ bonds and 
one NNN J2 bond. The Shastry-Sunderland model finds 
a good experimental realization in the magnetic material 
SrCu(B03)2- A second such model, in which half of the 
J 2 bonds of the J1-J2 model are removed in an arrange- 
ment that leaves each lattice site connected by four NN 
Ji bonds and two NNN J 2 bonds, is the anisotropic trian- 
gular lattice (or Ji-J 2 ) model [25], which is also believed 
to well describe the layered magnetic material CS2CUCI4. 
A third such model, in which again half of the J 2 bonds 
of the J\-J 2 model are removed, but now in an arrange- 
ment that leaves half of the sites eight-connected (namely, 



by four NN J\ bonds and four NNN J2 bonds) and the 
other half four-connected by J\ bonds alone, is the Union 
Jack model [11 15011511152] . whose study we continue here 
by extending the situation to where the spins have spin 
quantum number s > ^ . Although all of the models men- 
tioned above show antiferromagnetic Neel ordering for 
small values of J2, their phase diagrams for larger values 
of J 2 display a wide variety of behaviour, including, two- 
dimensional quantum "spirals" , valence-bond crystals/solids, 
and spin liquids [33] . Thus, in the absence of any definitive 
overarching theoretical argument, the best way to under- 
stand this class of NN/NNN models on the square lattice 
is to treat each one on a case-by-case basis. 

While the general trend is that as the spin quantum 
number s is increased the effects of quantum fluctuations 
reduce, one also knows that there can be significant devi- 
ations from it. A particularly well-known example is the 
since-confirmed prediction of Haldanc that integer-spin 
systems on the linear chain would have a nonzero exci- 
tation energy gap, whereas half-odd-integer spin systems 
would be gapless [33]. Indeed, such deviations from gen- 
eral trends provide one of the main reasons why quantum 
spin-lattice problems still maintain such an important role 
in the general investigation of quantum phase transitions. 

For the past few decades, a great deal of attention has 
also been devoted to magnetic materials with spin-1 ions 
(see, e.g., Refs. [5S] and references cited therein). In this 
context we note the recent discovery of superconductivity 
with a transition temperature at T c « 26 K in the layered 
iron-based compound LaOFcAs, when doped by partial 
substitution of the oxygen atoms by fluorine atoms [36], 
Lapi-zF^FeAs, with x w 0.05-0.11. This has been fol- 
lowed by the rapid discovery of superconductivity at even 
higher values of T c (> 50 K) in a broad class of similar 
doped quaternary oxypnictide compounds. Enormous in- 
terest has thereby been engendered in this class of mate- 
rials. Of particular relevance to the present work are the 
very recent first-principles calculations |37j showing that 
the undoped parent precursor material LaOFeAs is well 
described by the spin-1 J\-J 2 model on the square lattice 
with Ji > 0, J 2 > 0, and J2/J1 ~ 2. Broadly similar con- 
clusions have also been reached by other authors [38]. It 
is clear, therefore, that the theoretical study of 2D quan- 
tum magnets with s = 1 (or s > 1) are worthy of pursuit 
in their own right, as well as to highlight differences with 
their s = \ counterpart. 

In this article we continue the study, begun in Ref. [1] 
for the s = \ case, of magnetic ordering in the Union Jack 
model, by now extending the discussion to the two cases 
where the spins have either s = 1 or s = |. The model 
itself is discussed in more detail in Sec. [2J below. For the 
s = h case the model has been studied previously using 
SWT [501151] and SE techniques [52], as well as by the 
CCM P used here. However, to our knowledge, no results 
have yet been reported for the cases with s > ^. 

As in the case of the spin-half anisotropic triangular 
lattice (or Ji-J 2 ) model, it was shown [3M5T1I32] that NN 
Neel order for the s = | Union Jack model persists until a 
critical value of the frustrating NNN ( J2) bonds. However, 
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Fig. 1: Union Jack model; — J\\ J2. (a) Canted state; 

(b) semi-striped state. The unit cell is a square of side 
length 2. 



in contrast to the case of the anisotropic triangular lattice 
model, there exists a ferrimagnctic ground state in which 
spins on the eight-connected sites cant at a nonzero angle 
with respect to their directions in the corresponding Neel 
state. This model thus exhibits an overall magnetic mo- 
ment in this regime, which is quite unusual for spin-^ 2D 
materials with only Heisenberg bonds and which therefore 
preserve (spin) rotational symmetries in the Hamiltonian. 
This model also presents us with a difficult computational 
task in order to simulate its properties. Here we wish to 
continue to study this model for the s = 1 and s = 3/2 
cases using the CCM, which has consistently been shown 
to yield insight into a wide range of problems in quantum 
magnetism. 



other. Henceforth we set J\ = 1. We consider the model 

equivalently defined by the Union Jack geometry in which 

there are two sorts of sites, namely the A sites with eight 

NN sites and the B sites with four NN sites, as shown in 
Fig. Ha). 

Considered as a classical model, (and thus correspond- 
ing to the quantum case in the limit where the spin quan- 
tum number s — > 00), the Union Jack model has only 
two ground-state (gs) phases as the parameter k is var- 
ied over the range (0, 00). A simple variational analysis 
for the classical model reveals that for < k < \ the gs 
phase is Neel-ordered, exactly as for the full J1-J2 model. 
The Neel ordering induced by the J\ bonds acting alone 
is hence preserved as the strength of the competing J2 
bonds is increased, until the critical value k^: 1 = 0.5 is 
reached. For k > k^} a new phase of lower energy emerges, 
just as in the full J1-J2 model. For the full J1-J2 model 
that new phase is a classical striped state in which alter- 
nate rows (or columns) of spins are arranged antiparallel 
to one another, whereas the new classical gs phase for the 
Union Jack model is the canted ferrimagnetic state shown 
in Fig. [IJa) . In the canted state the spins on each of the 
alternating Ai and A2 sites of the A-sublattice arc canted 
respectively at angles (tt =p 0) with respect to those on the 
B sublattice, all of the latter of which point in the same 
direction. On the A sublattice each site Ai has four NN 
sites A2 , and vice versa. The angle between the NN spins 
on the A sublattice is thus 2cf>. 

The energy of the above canted state for the classical 
model is thus 



Ns 2 (k cos 20 - 2 cos < 



(2) 



where Ji = 1 and N — > 00 is the number of sites. The 
energy is thus extremized when 



2 The model 

The Hamiltonian of the Union Jack model considered here 
is given by 



H = Ji 



(id) 



J 2 ]>> 

[i,k] 



Sfc 



(1) 



where the operators s, = (sf , s^sf) are the quantum spin 
operators on lattice site i with sf = s(s + 1) where, for 
the cases considered here, s = On the underlying 

2D square lattice the sum over runs over all distinct 
NN bonds with strength J\, while the sum over [i, k] runs 
over only half of the distinct NNN diagonal bonds hav- 
ing strength J2 and with only one diagonal bond on each 
square plaquette as arranged in the pattern shown explic- 
itly in Fig. [TJ The unit cell is thus the 2x2 square shown 
in Fig. [TJ (We note that, by contrast, the J1-J2 model 
discussed above includes all of the diagonal NNN bonds 
on the square lattice, and its unit cell is thus the lxl 
square.) We again consider the case where both sorts of 
bonds are antiferromagnetic, J\ > and J2 = kJi > 0, 
and thus act to compete against (or to frustrate) each 



sm < 



[1 - 2k cos 6) = . 



(3) 



When k < K,f = 0.5, the lowest energy corresponds to 
sin0 = and hence to the Neel state. By contrast, when 
0.5 the lowest energy solution is the canted 



k > K 



:1 — 



state with 



Thus the classical gs energy is given by 



E c 



Ns 2 (k - 2) ; 



K< K C c l = 0.5 . 

k> K c } = 0.5 . 



(4) 



(5) 



The classical phase transition at k = n c c = 0.5 is of contin- 
uous (second-order) type with the gs energy and its first 
derivative with respect to n both continuous functions of 
k, although there are finite discontinuities in the second- 
and higher-order derivatives at k = k^}. 

The total magnetization per site in the canted phase 
for the classical model is m cl = |s[l — (2k) and the 
model thus exhibits ferrimagnctism in this phase. Whereas 
fcrrimagnetism more commonly occurs when the individ- 
ual ionic spins have different magnitudes on different sub- 
lattices, it arises here in a case where when the spins all 
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have the same magnitude and all the interactions are anti- 
ferromagnetic in nature, but the frustration between them 
acts to produce an overall magnetization. The total mag- 
netization to vanishes linearly asc^ from the canted 
phase and then remains zero in the Neel phase for k < nf. 
The spontaneous breaking of the spin rotation symmetry 
is also reflected by the vanishing of the energy gap on 
both sides of the transition. Clearly on both sides of the 
transition the translation symmetry of the lattice is also 
broken. 

As previously for the s = \ case [I] , our aim here is to 
give a fully microscopic analysis of the Union Jack model 
for the quantum case where the spins all have spin quan- 
tum number either s = 1 or s = |. Our goal is to map 
out the zero-temperature (T = 0) phase diagrams for both 
cases, including the positions and orders of any quantum 
phase transitions that emerge. In particular we investigate 
the quantum analogues of the classical Neel and canted 
phases and calculate the effect of quantum fluctuations on 
the position and nature of the transition between them. 
We also aim to investigate, for particular regions of the 
control parameter k, whether the quantum fluctuations 
may favour other phases, which have no classical counter- 
parts. One such possible candidate is discussed below. 

For the classical (s — > oo) model the k — > oo limit 
corresponds to a canting angle <fi — > ^ir, such that the 
spins on the A sublattice become Neel-ordered, as is ex- 
pected. The spins on the antifcrromagnctically-ordcrcd 
A sublattice are thus orientated at 90° to those on the 
ferromagnetically-ordered B sublattice in this limit. In re- 
ality, of course, there is complete degeneracy at the clas- 
sical level in this limit between all states for which the 
relative ordering directions for spins on the A and B sub- 
lattices are arbitrary. Clearly the exact ft — > oo limit of 
the Union Jack model with spins having any value of the 
spin quantum number s should also comprise decoupled 
antifcrromagnetic and ferromagnetic sublattices. However, 
one might now expect that this degeneracy in the relative 
spin orientations between the two sublattices is lifted by 
quantum fluctuations by the well-known phenomenon of 
order by disorder |39| . Just such a phase is known to exist 
in the full spin-i J\-Ji model for values of J2/J1 > 0.6, 
where it is the so-called collinear striped phase in which, 
on the square lattice, spins along (say) the rows in Fig. [T] 
order ferromagnetically while spins along the columns and 
diagonals order antifcrromagnetically. We have also shown 
how such a striped state is stabilized by quantum fluctu- 
ations for values of J^/Ji ^ 1-8 for the spin--| anisotropic 
triangular lattice model (or J\~J' 2 model) [29] . 

The existence of the striped state as a stable phase for 
large values of the frustration parameter for both the spin- 
\ J1—J2 and J\-J'i models above is a reflection of the well- 
known fact that quantum fluctuations favour collinear or- 
dering. In both cases the order-by-disorder mechanism 
favours the collinear state from the otherwise infinitely 
degenerate set of available states at the classical level. For 
the present Union Jack model the corresponding collinear 
state that might perhaps be favoured by the order by dis- 
order mechanism is the so-called semi-striped state shown 



in Fig. db) where the A sublattice is now Neel-ordered 
in the same direction as the B sublattice is ferromagnet- 
ically ordered. Alternate rows (or columns) are thus fer- 
romagnetically and antifcrromagnetically ordered in the 
same direction. We investigate the possibility below that 
if such a semi-stripe-ordered phase may be stabilized by 
quantum fluctuations at larger values of ft for either of the 
cases s = 1 or s = |, in order to compare with the earlier 
s = I case pQ. 



We note that for the s = ^ case our own CCM calcu- 
lations [T] provided strong evidence that the canted fcr- 
rimagnetic phase becomes unstable at large values of the 
frustration parameter k. In view of that observation we 
also used the CCM for the s = \ case with the collinear 
semi-stripe-ordered ferrimagnetic state as a model state. 
We found tentative evidence, based on the relative ener- 
gies of the two states, for a second zero-temperature phase 
transition between the canted and semi-stripe-ordered fer- 
rimagnetic states at a larger critical value of n C2 « 125±5, 
as well as firm evidence for a first phase transition between 
the Neel antifcrromagnetic phase and the canted ferrimag- 
netic phase at a critical coupling k Ci = 0.66 ± 0.02. Our 
prediction for k C2 , however, was based on an extrapolation 
of the CCM results for the canted state into regimes where 
the solutions have already become unstable and the CCM 
equations based on the canted state as model state have no 
solutions at any level of (LSUBn) approximation beyond 
the lowest (with n = 2). The prediction for k C2 was thus 
less reliable than that for k Ci , although our results showed 
clear evidence that, if the second transition at k = n C2 
does exist, it should be of first-order type. By contrast, 
the transition at k = k Ci for the s = ^ case was found to 
be an interesting one. As in the classical (s — > 00) case, 
the energy and its first derivative were seen to be contin- 
uous (within the errors inherent in our approximations), 
thus providing a typical scenario of a second-order phase 
transition, although a weakly first-order one could not be 
excluded. Nevertheless, the average on-site magnetization 
was seen to approach a nonzero value M Cl = 0.195 ±0.005 
on both sides of the transition, which is more typical of 
a first-order transition. The slope, dM/dn, of the order 
parameter curve as a function of the coupling parameter 
«, also appeared to be continuous, or very nearly so, at 
the critical point n Cl . Thus, all of the evidence shows that 
for the s = \ case the transition between the Neel and 
canted phases is a subtle one. A particular interest here is 
to compare and contrast the corresponding transition(s) 
between the s = \ and the s > \ models. 



We first briefly describe the main elements of the CCM 
below in Sec. [21 where we also discuss the approximation 
schemes used in practice for the s = \ case and the s > | 
cases. Then in SecHJwe present our CCM results based on 
using the Neel, canted and semi-striped states discussed 
above as model states (or starting states). We conclude in 
Sec. [S] with a discussion of the results. 
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3 The coupled cluster method 

It is widely recognized nowadays that the CCM (see, e.g., 
Rcfs. jHElH] an d references cited therein) employed here is 
one of the most powerful and most versatile modern tech- 
niques in quantum many-body theory. It has been suc- 
cessfully applied to various quantum magnets (see, e.g., 
Refs. [^[6l[71pll[T5l[T6l[29l[¥)l[4Tl[42] . and is particularly ap- 
propriate for studying frustrated systems, for which the 
main alternative methods are often only of limited use- 
fulness. For example, QMC techniques are particularly 
plagued by the sign problem for such systems, and the 
ED method is often restricted in practice, particularly for 
s > 1/2, to such small lattices that it is often insensitive 
to the details of any subtle phase order present. 

We now briefly describe the CCM means to solve the 
ground-state (gs) Schrodingcr ket and bra equations, H \>P) 
E\W) and (§\H = respectively (and see Refs. 011 

I4l40ll41ll4"2"] for further details). The first step in imple- 
menting the CCM is always to choose a normalized start- 
ing state or model state \<P) on top of which to incor- 
porate later in a systematic fashion the multi-spin cor- 
relations contained in the exact ground states I'P) and 
(<?|. More specifically, the CCM employs the exponen- 
tial parametrizations, \W) = e s |<£) and (i^l = (<P\Se~ s . 
These states arc chosen with a normalization such that 
= 1 [i.e., with (#| = and with \W) 

itself satisfying the intermediate normalization condition 
(<P\ X I / ) = 1 = ($\<1>)- The correlation operator S is ex- 
pressed as S — ^2i^ SjC^ and its counterpart is S — 
! + E/^o SiCJ. The operators Cj = (Cj)\ with C+ = 1, 
have the property that (<£|C/ = = CJ\$); V/ ^ 0. They 
form a complete set of multi-spin creation operators with 
respect to the model state \<P). The index / is a set-index 
that stands for the set of lattice sites whose spin projec- 
tions on the quantization axis are changed with respect to 
their values in the model state The ket- and bra-state 
correlation coefficients (Si, Si) arc calculated by requir- 
ing the gs energy expectation value H = (\P\H\ils) to be a 
minimum with respect to each of them. This immediately 
yields the coupled set of equations (&\Cj c~ s Hc s \<P) = 
and ($\S(e~ s Hc s -E)Cf\$) = ; V/ ^ 0, which we solve 
in practice for the correlation coefficients (Si, Si) within 
specific truncation schemes described below, by making 
use of parallel computing routines [43) . 

It is important to note that while the above CCM 
parametrizations of \W) and (<?| are not manifestly Hermi- 
tian conjugate, they do preserve the important Hellmann- 
Feynman theorem at a//levels of approximation (viz., when 
the complete set of many-particle configurations {/} is 
truncated [3]. Furthermore, the amplitudes (Si, Si) form 
canonically conjugate pairs in a time-dependent version of 
the CCM, in contrast with the pairs (Si,Sj) that come 
from a manifestly Hermitian conjugate representation for 
(W\ = ((^|e sl e s |<£)) _1 (^|e sl , which are not canonically 
conjugate to each other [3]. 

In order to treat each lattice site on an equal footing 
we perform a mathematical rotation of the local spin axes 



on each lattice site, such that every spin of the model state 
aligns along its negative z-axis. Henceforth our description 
of the spins is given wholly in terms of these locally de- 
fined spin coordinate frames. In particular, the multi-spin 
creation operators may be written as Cf = sf sf ■ ■ ■ sf , 
in terms of the locally defined spin-raising operators sf = 
sf + sf on lattice sites i. Having solved for the multi-spin 
cluster correlation coefficients (Si, Si) as described above, 
we may then calculate the gs energy E from the relation 
E = (<P\c~ s Hc s \<P), and the average gs on-site magnetiza- 
tion M from the relation M = J2?=i s iW) which 
holds in the rotated local spin coordinates. The quantity 
M is thus the magnetic order parameter, and it is just the 
usual sublattice (or staggered) magnetization for the case 
of the Neel state as CCM model state, for example. 

Although the CCM formalism is clearly exact if a com- 
plete set of multi-spin configurations {/} with respect to 
the model state \<1>) is included in the calculation of the 
correlation operators 5" and S, in practice it is necessary to 
use systematic approximation schemes to truncate them 
to some finite subset. In our earlier paper on the s = \ ver- 
sion of the present model [T] , we employed, as in our previ- 
ous work [2113:1,42,44,43], the localized LSUB?i scheme in 
which all possible multi-spin-ffip correlations over differ- 
ent locales on the lattice defined by n or fewer contiguous 
lattice sites are retained. 

However, we note that the number of fundamental 
LUBn configurations for s = 1 and s = | becomes ap- 
preciably higher than for s = ^, since each spin on each 
site i can now be flipped up to 2s times by the spin-raising 
operator sf . Thus, for the cases of s = 1 and s = | it is 
more practical, but equally systematic, to use the alterna- 
tive SUBn-m scheme, in which all correlations involving 
up to n spin flips spanning a range of no more than m ad- 
jacent lattice sites are retained [4"ll4"0]. We then set m = n, 
and hence employ the so-called SUBn-n scheme. More 
generally, the LSUBm scheme is thus equivalent to the 
SUBn-m scheme for n = 2sm for particles of spin s. For 
s = i, LSUBn = SUBn-n; whereas for s = 1, LSUBn = 
SUB2n-n, and for s = §, LSUBn = SUB3n-n. The num- 
bers of such fundamental configurations (viz., those that 
are distinct under the symmetries of the Hamiltonian and 
of the model state |^}) that are retained for the Neel and 
semi-striped states of the current s = 1 and s = | models 
at various SUBn-n levels are shown in Table [TJ We note 
that the distinct configurations given in Table [T] are de- 
fined with respect to the Union Jack geometry described 
in Sec. [2j in which the B sublattice sites of Fig. [lja) are 
defined to have four NN sites, and the A sublattice sites 
are defined to have the eight NN sites joined to them by 
either J\ bonds or Ji bonds. If we chose instead to work in 
the square-lattice geometry, by contrast, each site would 
have four NN sites. 

Although we never need to perform any finite-size scal- 
ing, since all CCM approximations are automatically per- 
formed from the outset in the TV — > oc limit, we do need 
as a last step to extrapolate to the n —y oo limit in the 
truncation index n. We use the same well-tested scaling 
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Table 1: Number of fundamental SXJBn-n configurations 
(Nf) for the semi-striped and canted states of the spin-1 
and spin-| Union Jack models, based on the Union Jack 
geometry defined in the text. 



s = 1 s = | 

Method Nf Nf 

semi-striped canted semi-striped canted 

SUB2-2 3 9 3 9 

SUB4-4 115 556 115 618 

SUB6-6 7826 52650 9862 68365 



the two minima have equal depth, plus one or two insta- 
bility points ft^ and Ki 2 where one or other of the minima 
(at 4> = and <p ^ respectively) disappears. By contrast, 
a second-order phase transition might manifest itself via 
a one-minimum structure for the gs energy as a function 
of </>, in which the single minimum moves smoothly and 
continuously from the value </> = for all values of k < k Ci 
to nonzero value 7^ for k > k Ci . 

Results for the gs energy per spin calculated in the 
SUB4-4 approximation based on the canted state as the 
CCM model state are shown in Fig. [3] as a function of the 
canting angle. Curves such as those shown in Fig. [5] for 
the SUB4-4 case show that what happens for this model 
at this level of approximation is that for k < 0.592 for 
s = 1 and for k < 0.560 for s = |, the only minimum in 
the gs energy is at (j> = (Neel order). As these values are 
approached from below the SUB4-4 energy curves become 
extremely flat near = 0, indicating the disappearance at 
= of the second derivative d 2 E/dcf> 2 (and possibly also 
of one or more of the higher derivatives d n E/d(f> n with 
n > 3), as well as of the first derivative dE/dcj). Then, for 
all values K > 0.592 for s = 1 and k > 0.560 for s = |, 
the SUB4-4 curves develop a minimum at a value 7^ 
which is also the global minimum. The state for 7^ is 
thus the quantum analogue of the classical canted phase. 
The fact that the antifcrromagnetic Neel order survives 
into the classically unstable regime is another example of 
the well-known phenomenon that quantum fluctuations 
tend to promote collincar order in magnetic spin-lattice 
systems, as has been observed in many other such cases 
(see e.g., Ref. [4"T1I46|). Thus, this collinear Neel-ordercd 
state survives into a region where classically it becomes 
unstable with respect to the non-collinear canted state. 
As expected, as s is increased the value of k Ci decreases 
towards the classical value 0.5 

A detailed inspection of the curves shown in Fig. [3] for 
various SUBn-?i approximation shows that the crossover 
from one minimum (0 = 0, Neel) solution to the other 
(0 7^ 0, canted) appears to be continuous for both the 
s = 1 and s = | cases indicating a second-order transi- 
tion according to the above scenario. Thus, based on the 
evidence presented so far of the gs energies of the Neel 
and canted phases, it would appear that the transition at 
k = k Ci between these two phases is second-order. 

Table [2] shows the critical values k^ubh-h & ^ which the 
transition between the Neel and canted phases occurs in 
the various SUBn-n approximations shown in Fig. [31 In 
the past we have found that a simple linear extrapolation, 
K SUBn-?i = cto + ain^ 1 , yields a good fit to such criti- 
cal points. This seems to be the case here too, just as for 
the spin-i case pQ. The corresponding "SUBoo" estimates 
from the SUBn-n data in Tableware K C1 = 0.585 ± 0.002 
for s = 1, and k Ci = 0.557 ± 0.001 for s = |, where the 
quoted errors are simply a combination of the standard 
deviations from the fits and of the computational uncer- 
tainties associated with the K suBn-n p 0m ^ s themselves. 
We also present other independent estimates of n Cl be- 
low. 



laws as for the s = 5 model for the gs energy per spin 
E/N and the average gs on-site magnetization M, 

E/N = a + a x nr 2 + a 2 ?i~ 4 , (6) 

M = bn + hn- 1 + b 2 nr 2 . (7) 



4 Results 

We present results of CCM calculations for both the spin- 
1 and spin- 1 Union Jack models with the Hamiltonian 
of Eq. ([T]), for given parameters (Ji = 1, J2), based re- 
spectively on the Neel, canted and semi-striped states as 
CCM model states. Our computational power is such that 
we can perform SUBn-77 calculations for each model state 
with n < 6. 



4.1 Neel state versus the canted state 

Results are first presented that are obtained using the 
Neel and canted model states. While classically we have 
a second-order phase transition from Neel order (for n < 
k^}) to canted order (for k > nf), where k = J2/J1, at a 
value tif = 0.5, using the CCM we find strong indications 
of a shift of this critical point to a higher value k Ci 0.58 
in the spin-1 case and k Ci « 0.545 in the spin-| quantum 
case as we explain in detail below. These may be compared 
with the even higher value k Ci ks 0.66 found previously [T] 
for the spin-i case. 

Thus, just as for the spin-^ case, curves such as those 
in Fig. [2] show that the Neel model state (<fi = 0) gives 
the minimum gs energy for all values of n < k Ci where 
k Ci = K SUBn-n j s a i s0 dependent on the level of SUBn- 
n approximation, as we see clearly in Fig. [3J Conversely, 
for k > k Ci the minimum in the energy is found to oc- 
cur at a value </> 7^ 0. If we consider the canting angle 
<t> itself as an order parameter (i.e., <fi — for Neel or- 
der and (j) 7^ for canted order) a typical scenario for a 
first-order phase transition would be the appearance of a 
two-minimum structure for the gs energy as a function of 
<f>. If we therefore admit such a scenario, in the typical case 
one would expect various special points in the transition 
region, namely the phase transition point k c , itself where 
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Fig. 2: Ground-state energy per spin of the spin-1 and spin-| Union Jack Hamiltonian of Eq. ([!} with J\ = 1. using 
the SUB4-4 approximation of the CCM with the canted model state, versus the canting angle <f>, for various values 
of J2. For the case of s = 1, for J2 < 0.592 in this approximation the minimum is at <fr = (Neel order), whereas 
for J2 > 0.592 the minimum occurs at <f> = 0suB4-4 7^ 0, indicating a phase transition at J2 ~ 0.592 in this SUB4-4 
approximation. For the case of s = |, for J2 < 0.560 in this approximation the minimum is at <f) = (Neel order), 
whereas for J2 > 0.560 the minimum occurs at <fi = </>suB4-4 7^ 0, indicating a phase transition at J2 f=s 0.560 in this 
SUB4-4 approximation. 



Figure [2] also displays the feature that for certain val- 
ues of J2 with Ji = 1 (or, equivalently, k) CCM solutions 
at a given SUBn-n level of approximation (viz., SUB4-4 

Table 2: The critical value K SUBn-n a ^ wn j cn transi- 
tion between the Neel phase (</> = 0) and the canted phase 
(4> 7^ 0) occurs in the SUBn-n approximation using the 
CCM with (Neel or) canted state as model state. 



Method 


s = 1 




suBn-n 


„SUB?i-n 


SUB2-2 


0.598 


0.563 


SUB4-4 


0.592 


0.560 


SUB6-6 


0.589 


0.559 


SUBoo 


0.585 ± 0.002 


0.557 ±0.001 



in Fig. [2]) exist only for certain ranges of the canting angle 
(j>. For example, for the pure square-lattice HAF (k = 0) 
the CCM SUB4-4 solution based on a canted model state 
only exists for < < 0.1827T when s = 1, and for 
< 4> < 0.1667T when s = |. In this case, where the 
Neel solution is the stable ground state, if we attempt to 
move too far away from Neel collinearity the CCM equa- 
tions themselves become "unstable" and simply do not 
have a real solution. Similarly, we see from Fig. [2] that 
for k = 1.5, for example, the CCM SUB4-4 solution ex- 
ists only for 0.2947T < </> < 0.4467T when s = 1 and for 
0.302-7T < (f> < 0.4467T when s = |. In this case the stable 
ground state is a canted phase, and now if we attempt 
either to move too close to Neel collinearity or to increase 
the canting angle too close to its asymptotic value of 7r/2, 
the real solution terminates. 

Terminations of CCM solutions like those discussed 
above are very common and have been very well doc- 
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(a) s = 1 (b) s = 3/2 

Fig. 3: The angle (j>svBn-n that minimizes the energy -EsuBn-n(0) of the (a) spin-1 and (b) spin-| Union Jack 
Hamiltonian of Eq. ([T]) with J\ = 1, in the SUBn-n approximations with n = {2,4,6}, using the canted model state, 
versus J2 . The corresponding classical result (f> c \ from Eq. Q is shown for comparison. Wc find in the SUBrt-n quantum 
case with n > 2 a second-order phase transition (e.g., for SUB6-6 at J 2 ~ 0.589 for s = 1 and J2 ~ 0.559 for s = |). 
In the classical case there is a second-order phase transition at J2 = 0.5. 



umented [4]. Such termination points always arise due 
to the solution of the CCM equations becoming complex 
there. Beyond such points there exist two branches of en- 
tirely unphysical complex conjugate solutions [4]. In the 
region where the solution reflecting the true physical so- 
lution is real there actually also exists another (unstable) 
real solution. However, only the shown branch of these two 
solutions reflects the true (stable) physical ground state, 
while the other branch does not. The physical branch is 
usually easily identified in practice as the one which be- 
comes exact in some known (e.g., perturbative) limit. It 
then meets the corresponding unphysical branch at some 
termination point (with infinite slope on Fig. [5]), beyond 
which no real solutions exist. The LSUB?i or SUBn-n 
termination points are themselves also reflections of the 
quantum phase transitions in the real system, and may be 
used to estimate the position of the phase boundary [3], 
although we do not do so for this first critical point since 
we have more accurate criteria, one of which has already 
been discussed above. Another will be discussed below. 

Before doing so, however, we wish to give some further 
indication of the accuracy of our results. Thus in Table [3] 
we show data for the cases of the spin-1 and spin-| HAF 
on the square lattice (corresponding to the case k = of 
the present Union Jack model). We present our CCM re- 
sults in various SUBn-n approximations (with 2 < n < 6) 
based on the Union Jack geometry using the Neel model 
state. Results are given for the gs energy per spin E/N, 
and the magnetic order parameter M, which in this case 
is simply the staggered magnetization. We also display 
our extrapolated (n — > 00) results using the schemes of 
Eqs. © and with the data set n = {2,4,6}. The 
results are robust, and for comparison purposes wc also 
show for the spin-1 case, the corresponding results us- 
ing a spin wave theory (SWT) technique |47j and from 
a linked-cluster series expansion (SE) method [15]. Our 



Table 3: Ground-state energy per spin E/N and magnetic 
order parameter M (i.e., the average on-site magnetiza- 
tion) for the spin-1 and spin-| square-lattice HAF. We 
show CCM results obtained for the Union Jack model with 
J\ = 1 and J2 = using the Neel model state in various 
CCM SUBn-n approximations defined on the Union Jack 
geometry described in Sec. [2] We compare our extrapo- 
lated (n — > 00) results using Eqs. © and (J7]) and our 
SUBn-n data sets with other calculations. 



Method 



E/N 



M 



E/N 



M 



SUB2-2 
SUB4-4 
SUB6-6 
SUBoo 
SWT a 
SE 6 



-2.29504 
-2.32033 
-2.32537 
-2.3295 
-2.3282 
-2.3279(2) 



0.9100 
0.8687 
0.8488 
0.800 
0.8043 
0.8039(4) 



-4.94393 1.4043 

-4.97758 1.3618 

-4.98352 1.3423 

-4.9882 1.295 



a SWT (Spin Wave Theory) for square lattice g7] 
b SE (Series Expansion) for square lattice [48] 



own extrapolated results are in good agreement with these 
results and own previous CCM results [4"9"| . 

The CCM results for the gs energy per spin are shown 
in Fig. 2] for various SUBn-n approximations based on 
the canted (and Neel) model states, with the canting an- 
gle ^suBn-n chosen to minimize the energy SsuBn-n(^), 
as shown in Fig. [31 We also show separately the extrap- 
olated (SUBoo) results obtained from Eq. © using the 
data set n = {2, 4, 6} shown. As is expected from our pre- 
vious discussion the energy curves themselves show very 
little evidence of the phase transition at K = n Cl , with 
the energy and its first derivative seemingly continuous 
for both the spin-1 and spin-| cases. 
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(a) s = 1 (b) s = 3/2 

Fig. 4: Ground-state energy per spin versus Ji for the Neel and canted phases of the (a) spin-1 and (b) spin-| Union 
Jack Hamiltonian of Eq. (0 with J\ = 1. The CCM results using the canted model state are shown for various 
SUB?i — ?i approximations with n = {2, 4, 6} with the canting angle </> = </>LSUBn that minimizes -EsuBn-n (</>)• We also 
show the n — > oo extrapolated result from using Eq. ©. 
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Fig. 5: Ground-state magnetic order parameter (i.e., the average on-site magnetization) versus Ji for the Neel and 
canted phases of the (a) spin-1 and (b) spin-| Union Jack Hamiltonian of Eq. (0 with J\ = 1. The CCM results 
using the canted model state are shown for various SUBn-n approximations with n = {2, 4, 6}, with the canting angle 
4> = 0suBn-n that minimizes £'suBn-n(0)- We also show the n — > oo extrapolated result from using Eq. 0. 



The transition between the Neel and canted phases 
is observed much more clearly in our corresponding re- 
sults for the gs magnetic order parameter M (the aver- 
age on-site magnetization) shown in Fig. [SJ For the raw 
SUB?i-n data we display the results for the Neel phase 
only for values of re < K SUBn-n f or c i ar ity. However, the 
extrapolated (SUBoo) results for the Neel phase are shown 
for a wider range of values of re, using the extrapola- 
tion scheme of Eq. (0 and the SUBn-n results based on 
the Neel model state. For the canted phase (for which 

r'SUBn— n 

i=- 0) we can clearly only show the extrapo- 
lated (SUBoo) results using Eq. (0, for regions of k for 
which we have data for all of the set n = {2,4,6}. We 
see from Table [5] that we are limited (by the SUB2-2 re- 
sults) to values re > k^ UB2_2 sa 0.598 for the s = 1 case, 
and k > k^ UB2_2 ps 0.563 for the s = | case. From 



the data shown in Fig. EJa) for the case of s — 1, sim- 
ple extrapolations of the SUBoo curve to lower values of 
re < k^ UB2_2 ps 0.598 using simple spline fits in re give a 
corresponding estimate of re Cl ps 0.580±0.015 at which the 
Neel and canted phases meet. The SUBn-n extrapolations 
yield a nonzero value for the average on-site magnetiza- 
tion of M ps 0.707 ± 0.003 at the phase transition point 
re Cl . Similarly for the case of s = | in Fig. [3Jb) , sim- 
ple extrapolations of the SUBoo curve to lower values of 
K < k^- TB2 ~ 2 ps 0.563 using simple spline fits in re give a 
corresponding estimate of re Cl ps 0.535±0.005 at which the 
Neel and canted phases meet. The SUBn-n extrapolations 
yield a nonzero value for the average on-site magnetiza- 
tion of M w 1.192 ± 0.002 at the phase transition point 
re Cl . Thus the evidence from the behaviour of the order 
parameter is that the transition at K C1 is a first-order one 
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for both the cases s = 1 and s = |, in the sense that 
the order parameter does not go to zero at k Ci , although 
it is certainly continuous at this point. The extrapolated 
curves also provide every indication that the derivative of 
the order parameter as a function of n is also continuous 
(or very nearly so) at K = k Ci for both the cases s = 1 
and s = §. 

In Fig. [5] we also show the corresponding extrapolated 
(SUBoo) results for the average on-site magnetization as 
a function of J 2 (with J\ = 1), or hence equivalently as a 
function of k, for both the A sites (A/a) and the B sites 
(Mb) of the Union Jack lattice. We recall that, as shown 
in Fig.[TJ each of the A and B sites is connected to four NN 
sites on the square lattice by J\ bonds, whereas each of the 
A sites is additionally connected to four NNN sites on the 
square lattice by J 2 bonds. The extrapolations are shown 
in exactly the same regions, and for the same reasons, as 
those shown in Fig. [5j 

We also comment briefly on the large- J 2 behaviour of 
our results for the canted phase. (We note that for compu- 
tational purposes it is easier to re-scale the original Hamil- 
tonian of Eq. (p} by putting J 2 = 1 and considering small 
values of J\.) For the case of s = 5 which we studied 
previously [I], the most interesting feature of the CCM 
results using the canted state as model state is that in all 
LSUBn approximations with n > 2 a termination point 
reached, beyond which no real solution can be 
found. For example, we found [T] the values k^ SUB4 « 80 
and k^ SUB6 sa 80. This is a first indication that the canted 
state becomes unstable at very large values of k against 
the formation of another (as yet unknown) state, as we dis- 
cuss further in Sec. l4.2l bclow. By contrast, in the present 
cases of s = 1 and s = §, the canted state shows no 
sign of any instability for any SUBn-n approximation with 
n = {2, 4, 6} for all values of n < 100 studied. 

Simple extrapolations to the n — » 00 limit of the gs 
energy using the data at K < 100 show that at large 
J 2 values we have E/N — > — 1.1652 J2 for s — 1 and 
E/N — 2.4943J2 for s — I . These numerical coefficients 
are almost exactly half of the values quoted in Table [3] for 
the case J 2 = 0. This is as expected since both the k — > 
and the n — > 00 limits of the Union Jack model are the 
square-lattice HAF, where in the latter case the square 
lattice contains only half the original sites, namely the A 
sites. Similarly, the extrapolated SUBoo values at larger 
values of k = 100 for the on-site magnetization on the A 
sites are M A 0.807 for s = 1 and M A -> 1.303 for s = §. 
These values are again in very good agreement with those 
shown in Table [3] for the J 2 = limit. The corresponding 
asymptotic values for the B-site magnetization are consis- 
tent with Mb — > 1.0 and Mb — > 1.5, as expected for large 
values of J 2 . 

4.2 Canted state versus the semi-striped state 

For the Union Jack model considered here, but for the case 
s = 5 studied previously, we predicted [1] a second phase 
transition at k = k C2 ps 125 ±5 using the same CCM tech- 
nique as used here. At this upper transition the ferrimag- 



nctic quantum canted phase shown in Fig. QJa) gives way 
to the quantum ferrimagnetic semi-striped phase shown 
in Fig. [IJb), such that for value n > k C2 the semi-striped 
phase becomes lower in energy. In order to investigate the 
possible similar stabilization of the quantum semi-striped 
state at high values of the frustration parameter k, we have 
performed similar CCM calculations here, for the s = 1 
and s = I cases, based on the semi-striped state as model 
state. 

As already noted above, we find no sign of instability 
in the canted state, for either of the cases s = 1 or s = | , 
for any level of SUBn-n approximation with n < 6 at all 
values k < 100 investigated, unlike for the corresponding 
LSUBn approximations in the s = i case that terminated 
in this range for all values 2 < n < 6 investigated. Fur- 
thermore, we find no evidence at all that the semi-striped 
phase has lower energy than the canted phase for any value 
of k for cither of the case s = 1 or s = |. Nonetheless, 
simple extrapolations to the k — > 00 limit of the gs energy 
of the semi-striped phase, using the SUBn-n data with 
n = {2, 4, 6} at k < 1000, show that at large J2 values we 
have E/N -> -1.1646J 2 for s = 1 and E/N -> -2.4939J 2 
for s = I (with J 1 = 1). The corresponding asymptotic 
(k — > 00) values for the average gs on-site magnetization 
of the semi-striped state are Ma — >• 0.807 on the A sites 
and Mb — > 1.303 on the A sites and Mb — >• 1.5 on the 
B sites for the s = | case, just as for the corresponding 
large-K limits of the canted state, and as corresponds to 
the pure square-lattice HAF as already noted above. 



5 Discussion and conclusions 

In an earlier paper [T] we used the CCM to study the effect 
of quantum fluctuations on the zero-temperature gs phase 
diagram of a frustrated spin-^ Heisenberg antiferromagnet 
(HAF) defined on the 2D Union Jack lattice of Fig. Q] 
In the present paper we have extended the analysis to 
the two computationally more challenging cases where all 
the lattice spins have spin quantum number either s = 1 
or s = I. We have once again concentrated on the case 
where the NN Ji bonds are antiferromagnetic (Jj > 0) 
and the competing NNN J 2 = kJi bonds in the Union 
Jack array have a strength in the range < k < 00. 
On the underlying bipartite square lattice there are thus 
two types of sites, viz., the A sites that are connected 
to the four NN sites on the B sublattice with Ji bonds 
and to the four NNN sites on the A sublattice with J 2 
bonds, and the B sites that are connected only to the 
four NN sites on the A sublattice with Ji bonds. The 
K = limit of the model thus corresponds to the HAF on 
the original square lattice (of both A and B sites) , while 
the k — > 00 limit corresponds to the HAF on the square 
lattice comprised of only A sites. We have seen that at the 
classical level (corresponding to the case where the spin 
quantum number s — > 00) this Union Jack model has only 
two stable gs phases, one with Neel order for k < k^} = 0.5 
and another with canted ferrimagnetic order for k > nf. 
We have therefore first used these two classical states as 
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Fig. 6: Extrapolated curves (SUBoo) for the ground-state magnetic order parameters (i.e., the on-site magnetizations) 
Ma at sites A (joined by eight bonds to other sites) and Mb at sites B (joined by four bonds to other sites) of the 
Union Jack lattice [and see Fig. [Tfa)] versus J2 for the Neel and canted phases of the (a) spin-1 and (b) spin-| Union 
Jack Hamiltonian of Eq. ((TJ with J\ = 1. The CCM results using the canted model state are shown for various SUBn-n 
approximations (n = {2,4,6}) with the canting angle = 0suBn-n that minimizes fisuBn-n(^)- 



CCM model states to investigate the effects of quantum 
fluctuations on them. 

For the spin-1 model we find that the phase transition 
between the Neel antiferromagnetic phase and the canted 
ferrimagnetic phase occurs at the value K C1 = 0.580 ± 
0.015, whereas for the spin-| model we find that the phase 
transition occurs at k Ci = 0.545 ± 0.015. The evidence 
from our calculations is that the transition at 
subtle one. From the energies of the two phases it ap- 
pears that the transition is second-order, as in the clas- 
sical case. However, on neither side of the transition at 
k Ci does the order parameter M (i.e., the average on-site 
magnetization) go to zero. Instead, as k —> n Cl from ei- 
ther side, M -)■ 0.707 ± 0.003 for the s = 1 case and 
M — > 1.192 ± 0.002 for the s = | case, which are more in- 
dicative of a first-order transition. Furthermore, the slope 
dM/dn of the average on-site magnetization as a function 
of k also seems to be cither continuous or to have only a 
very weak discontinuity at K = k Ci in both of the cases 
s = 1 and s = |. 

In the case of the previously studied spin-i model [1] 
we found evidence for a second quantum phase transition 
at a value k = k C2 rj 125 ± 5 at which the quantum ferri- 
magnetic canted phase for k < k C2 yields to a lower-energy 
quantum ferrimagnetic semi-striped phase for k > n C2 . 
Our LSUBn results with n > 2 for the spin-i model gave 
clear evidence that the canted phase terminated at some 
value K^SUBn gucn tnat f or K > K LSUBn n0 rea i so l u tion 

based on the canted state as model state existed. These 
termination points were a preliminary signal of the actual 
phase transition at k C2 . In neither the s = 1 nor s = | 
cases considered here do we find any corresponding upper 
termination points for any SUBrt-n approximation based 
on the canted state as model state, for any value of n < 6 
and for any value of k up to the highest value k = 100 
investigated. Furthermore, for all approximation and for 
all values of k investigated the semi-striped phase always 



lies higher in energy than the quantum canted phase for 
both the s = 1 and s = | cases, unlike what was found [T] 
for the s = \ case. Thus, our results clearly show that 
only the s = i model has a second quantum phase tran- 
sition at k = k C2 between the two ferrimagnetic (canted 
and semi-striped) phases. No such transition is predicted 
to occur at a finite value of k for either of the s = 1 or 
s = § cases. 

To the best of knowledge no other studies of the Union 
Jack model for lattice spins with spin quantum number 
s > ^ have been performed, and hence we have no other 
results with which to compare. Nevertheless, the consis- 
tency of the present results and those of our previously 
studied spin-i case [T] give credence to our results. Fur- 
thermore, as has been noted elsewhere [7], high-order CCM 
results of the sort presented here have been seen to pro- 
vide very accurate and reliable results for a wide range 
of such highly frustrated spin-lattice models. Many pre- 
vious applications of the CCM to unfrustrated spin mod- 
els have given excellent quantitative agreement with other 
numerical methods (including exact diagonalization (ED) 
of small lattices, quantum Monte Carlo (QMC), and se- 
ries expansion (SE) techniques). A typical example is the 
HAF on the square lattice, which is the k = limit of the 
present model [T]. We have compared our own results for 
this k = case with both SWT and SE results, both for 
the spin-i case in our earlier work pQ and for the spin-1 
case here, and shown that there is excellent agreement in 
both cases. This adds further credence to the validity of 
our results at nonzero values of k. 
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